arXiv:1509.04596v2 [hep-ph] 25 Jan 2016 


JSNU/PHY-TH-2015 


Large non-factorizable contributions in 5 —> agao decays 

Defa Dou, Xin Liu^’’, and Jing-Wu Li‘^ 

School of Physics and Electronic Engineering, 

Jiangsu Normal University, Xuzhou, Jiangsu 221116, People’s Republic of China 

Zhen-Jun Xiao‘* 

Department of Physics and Institute of Theoretical Physics, 

Nanjing Normal University, Nanjing, Jiangsu 210023, People’s Republic of China 

(Dated: January 26, 2016) 

We investigate three tree-dominated B —>• aoao decays for the first time in the perturbative 
QCD(pQCD) approach at leading order in the standard model, with oq standing for the light 
scalar oo(980) state, which is assumed as a meson based on the model of conventional 
two-quark(gg) structure. All the topologies of the Feynman diagrams such as the non- 
factorizable spectator ones and the annihilation ones are calculated in the pQCD approach. 
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I. INTRODUCTION 


As we know, the nature of the light sealar states such as ao(980) is not yet well understood at 
both theoretical and experimental aspects. Also the identification and the classification of these 
light scalars remain as a long-standing puzzle(for latest review, see, e.g. [1]) to be resolved. How¬ 
ever, it is fortunate for the people that the light scalars as products in the heavy flavor meson de¬ 
cays have been detected, for example, D ^ SP, SV, B ^ SP, SV, even B ^ SS modes [1,2] 
with S, P, and V being the light scalar, pseudoscalar, and vector mesons, respectively, which 
will provide unique places and play very important roles on investigating the physical properties 
of light scalars. It is generally believed that the ongoing Large Hadron Collider(LHC) experi¬ 
ments can provide rich data on the B, Bg, and B^ meson decaying into light scalars. And more 
promisingly, the forthcoming Super-B/Belle-II factory scheduled in 2018 with a high luminosity 
> [3, 4] will produce much more events about the relevant decays. The studies on 

the above mentioned decays can also provide more constraints complementarily on the parameters 
in the standard model(SM), hint the exotic new physics beyond the SM, etc. 

In this work, we will investigate the CP-averaged branching ratios and the CP-violating asym¬ 
metries of the B —)■ ao(980)ao(980) decays by employing the perturbative QCD(pQCD) ap¬ 
proach [5-7] with the low energy effective Hamiltonian [8] in the SM. It should be noted that 
the ao(980) state here will be assumed as a meson in the model of conventional two-quark(q'g) 
structure. Moreover, hereafter, the ao(980) will be abbreviated as oq for the sake of simplicity 
throughout the paper. To our knowledge, heretofore, no other B ^ SS processes have been 
studied explicitly in the factorization approaches based on the QCD dynamics, apart from the 
Bu,d,s — ^ A'o (1430) J?o (1430) decays [9] by two of our authors (X. Liu and Z.J. Xiao). Because 
the scalar meson has either tiny or vanishing vector decay constant [10, 11], the contributions aris¬ 
ing from the factorizable emission diagrams in the B ^ SS decays are usually highly suppressed, 
which is dramatically different from the known B —)■ PP, PV, VV decays. In other words, for ex¬ 
ample, in contrast to the extensively investigated B ^ tttt decays, the large measured B —)■ aoao 
decay rates may indicate large non-factorizable spectator scattering and/or annihilation contribu¬ 
tions, which would hint some useful information on the i? —)■ tttt decays, the presently known puz¬ 
zle to be resolved, because they embrace the same components at the quark level. In the heavy B 
meson decays, the above mentioned large contributions from non-factorizable spectator and anni- 
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hilation diagrams are often eonsidered as the small ^ and/or negligible higher order or higher power 
eorreetions in the naive faetorization approaeh [13]. Therefore, the ehannels involving an emitted 
sealar state in the heavy flavor meson deeays are suggested to test the breaking effeets of the fae¬ 
torization assumption, e.g. [14]. Though the QCD improved faetorization approaeh [15, 16] going 
beyond the naive faetorization, the end-point singularities make it less predietive beeause the non- 
faetorizable speetator seattering eontributions and the annihilation ones have to be parametrized 
with the tunable parameters, whieh are always determined by the experimental measurements. As 
one of the popular faetorization approaehes based on the QCD dynamies, the pQCD approaeh 
involves no end-point singularities by retaining the parton transverse momentum Based on /cy 
faetorization theorem, the double logarithms arising from the overlap of soft and eollinear diver- 
genees generated in the radiative eorreetions are resummed into an important Sudakov faetor to 
suppress the long-distanee eontribution [17]. Armed with this pQCD approaeh, all the transition 
form faetors, the non-faetorizable speetator diagrams, and the annihilation diagrams are perturba- 
tively ealeulable, besides the faetorizable speetator diagrams. Note that, as far as the annihilation 
eontributions are eoneemed, soft-eollinear effeetive theory [18] and pQCD approaeh have an ex¬ 
tremely different effeet on the perturbative ealeulations [19, 20]. However, the predietions on the 
pure annihilation deeays based on the pQCD approaeh ean aeeommodate the experimental data 
well, for example, see Refs. [21-24]. We will therefore put the eontroversies aside and adopt the 
pQCD approaeh in our analyses. 

The paper is organized as follows. Seetion II is devoted to the analytie expressions for the 
deeay amplitudes of i? —?■ aoao modes in the pQCD approaeh. The numerieal results and phe- 
nomenologieal analyses on the CP-averaged branehing ratios and the CP-violating asymmetries of 
the eonsidered deeays are given in See. III. We summarize and eonelude in See. IV. 


II. PERTURBATIVE CALCULATIONS 


For the eonsidered B —)■ aoao deeays, the related weak effeetive Hamiltonian ifgff [8] can be 
written as 

G 

^ i=3 


HeS = 


V2 


10 


v:,v^i[CiWO';{p) + CMO'M] - v,iVulJ2Ci(ii)o,{p)] + h.c. , (i) 


* In fact, the cancelation of the decay amplitudes indeed occurred between the two non-factorizable spectator dia¬ 
grams in the B PP, PV, VV channels, for example, see Ref. [12]. 
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with the Fermi eonstant Gp = 1.16639 x 10“®GeV“^, the Cabibbo-Kobayashi-Maskawa(CKM) 
matrix elements V, and the Wilson eoeffieients Ci{jj) at the renormalization scale fi. The local 
four-quark operators Oi(z = 1, • • • , 10) are written as 

(1) current-current(tree) operators 

Oi = {daUp)v-A{upba)v-A , O 2 = {daUa)v-A{upbii)v-A 

(2) QCD penguin operators 


O 3 — {daba)v-A'^2^^Wl3)v-A , O 4 — , 

<?' <?' 

O5 = {daba)v-A'^2^^Wl3)v+A , Oq = ; 

q' q' 

(3) electroweak penguin operators 


( 3 ) 


Oj = -{daba)v-A'^eq>{q'0'^)v+A , Os = -{dJ)p)v-A'^eq'{q'0'^)v+A , 

9' (4) 

2 _ _^ 2 _ _^ 

O 9 = 2(daba)v-A'^eqi{q0'^)v-A , Olo = -{dabi3)v-A'^eq/{qpq'^)v-A ■ 

q' q' 

with the color indices a, fd and the notations {q'q')v±A = q'l^l{^ ± The index q' in the 

summation of the above operators runs through m, d, s, c, and b. The standard combinations Oj of 
Wilson coefficients are defined as follows, 

ai = C 2 + — , 02 = C*! + — , Oj = Gj H- -—{i = 3 — 10) . (5) 

where the upper(lower) sign applies, when i is odd(even). 

Similar to i? —)■ tttt decays [6, 25], there are eight types of diagrams contributing to i? —)■ aoOo 
modes at leading order(LO) in the pQCD approach, as illustrated in Fig. 1. They can be clas¬ 
sified into two types of topologies as emission and annihilation, respectively. And each kind of 
topology contains factorizable diagrams such as Fig. 1(a) and 1(b), in which a hard gluon con¬ 
nects the quarks in the same meson, and non-factorizable diagrams such as Fig. 1(c) and 1(d), in 
which a hard gluon attaches the quarks in two different mesons. By evaluating all these Feynman 
diagrams, one can obtain the decay amplitudes oi B ^ Uoao decays. Because the above men¬ 
tioned diagrams are the same as those in i? —)■ iTg ( 1430) -Aq (1430) modes [9], and also the light 
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FIG. 1. (Color online) Typical Feynman diagrams for B —)• oqOq decays at leading order in the pQCD 
approach. 


scalar mesons are considered, the formulas of B —)■ agOo decays are therefore same as those of 
B A'q (1430).^o( 1430) ones just by replaeing the wave funetions and input parameters eorre- 
spondingly. Hence the analytic formulas for the B —)■ agag decays are not explicitly presented in 
this paper. 

By taking various eontributions from the relevant Feynman diagrams into consideration, the 
total deeay amplitudes for three tree-dominated B —)■ aoao ehannels can then be read as, 

1. for B^ —)■ Oq Oq deeay mode. 


—y ®o'®o) — 


ClMnfs + C2Mnfa 


A, 


(C's + Cg)Mnfs + (^3 + 20^ — 


-Cio))Mnfa + (C 5 - + {2Cq + + (ag - -as) 


+ (oe + 


( 6 ) 


where A„ = V*^Vud and \t = V*^Vtd- We adopt F and M to denote the contributions from 
iy — A){V — A) operators in the faetorizable and non-factorizable diagrams, respeetively. 
Analogously, and are ehosen to denote the eontributions from (V — A) {V + A) 
operators, and and are taken to denote the eontributions from (S' — P){S + P) 
operators which result from the Fierz transformation of the (V — A){V + A) operators. 
The subseripts fs, nfs, fa, and nfa are the abbreviations for faetorizable emission, non- 
factorizable emission, faetorizable annihilation, and non-factorizable annihilation, respee- 
tively. 
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2. for 5+ —)■ Oq Oq decay mode, 


\/2vA(-B~*~ — y flo'flg) — 


(Cl + C2)Mnfs 


-Xf. 


1 , ^ ^ , ^P2 3 


+^CsM^}s 


;(Cr + + ;j(C9 + Cio)M„;, 


(V) 


3. for —)■ agOo decay mode, 


V2A{B^ a°a°) = A, 


C2{Mnfa — Mnfs) 




-(a, - \a»)Ff‘ + (Cs - ^(C, 


+ 3 Ci|)))Mn ^5 — + (C 3 + 2 C 4 — —(Cg — Ci(,j)Mnfa 

+(C, - + ( 2 Ce + + (a, - 


( 8 ) 


It is worth mentioning that the highly suppressed Ffs has been safely neglected in all of the above 
decay amplitudes for the considered B —)■ ogOg decays due to the either extremely small or van¬ 
ishing vector decay constant. Furthermore, based on the discussions of Ffa below Eq. (40) in 
Ref. [9], the factorizable annihilation contributions induced by the V ± A currents are therefore 
naturally absent because of the isospin symmetry between u and d quarks in the above analytical 
decay amplitudes. 


III. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will make theoretical predictions on the CP-averaged branching ratios and 
the CP-violating asymmetries for the B —)■ ogOg decay modes considered. In numerical calcula¬ 
tions, central values of the input parameters will be used implicitly unless otherwise stated. Firstly, 
we shall make several essential discussions on the input quantities. 


A. Input quantities 


For B meson, the distribution amplitude in the impact b space, with b being the conjugate space 
coordinate of transverse momentum kx, has been proposed [5-7], 


b) = Nbx^{ 1 — xY exp 


[ 1 

f XTUb a ^ 


2 

\ ) 

2 


(9) 
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where the normalization factor Nb is related to the decay constant fs through the following nor¬ 
malization condition, 


fn 

dx(j)B{x, b = 0) = 


2^/Wc 


( 10 ) 


with the color factor Nc = 3. The shape parameter Ub has been fixed at 0.40 GeV associated with 
Nb = 91.745 by using the rich experimental data on the B mesons with /b = 0.19 GeV based on 
lots of calculations of form factors and other well-known decay modes of B meson in the pQCD 
approach [5, 6 , 26]. 

For the light scalar oq, its leading twist light-cone distribution amplitude (paoix, /r) can be gen¬ 
erally expanded as the Gegenbauer polynomials [10, 27]: 


0ao(a;,/i) = “ ^)|/ao(/^) + /ao(/^) {2x - 1) 


( 11 ) 


where /ao(/^) and fao{lA^ and are the vector and scalar decay constants, Gegen¬ 

bauer moments, and Gegenbauer polynomials, respectively. For the vector and scalar decay con¬ 
stants, fao = fiaofao with fiao = = 0-985 GcV, whcrc rui and m 2 are the 

running current quark masses in the scalar oq. For neutral scalar Oq meson, which cannot be pro¬ 
duced by the vector current, the vector decay constant /„o = 0 is guaranteed by charge conjugation 
invariance. But the quantity fag = /ao/^ao remains finite. In fact, for the charged Og meson, the 
vector decay constant f^± also vanishes in the isospin limit. The reason is that f^± is proportional 
to the mass difference between the constituent d and u quarks, which will result in f^± being of or¬ 
der md — mu- Hence, the contribution from the first term in Eq. (11), namely, fag , can be neglected 
safely. In other words, the factorizable spectator diagrams could not contribute to 77 —?■ agao de¬ 
cays through the vector currents. We shall use the same light-cone distribution amplitudes for both 
neutral and charged oq mesons for simplicity in this paper. 

The values for scalar decay constant and Gegenbauer moments in the Og distribution amplitudes 
have been investigated at scale /i = 1 GeV [10]: 

/ao = 0.365 ± 0.020 GeV, =-0.93 ± 0.10 , ^3 = 0.14 ± 0.08 . (12) 


As for the twist-3 distribution amplitudes and cp'^^, we here adopt the asymptotic forms in 
our numerical calculations for simplicity [ 10 ]: 

“ 2x/2iVe'^““’ ~ 


ao 


(13) 
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The QCD scale (GeV), masses (GeV), and B meson lifetime(ps) are [1, 5, 6] 

= 0.250 , mw = 80.41 , ms = 5.2792 , = 4.8 ; 

tb+ = 1.643 , Tbo = 1.53 , rua, = 0.985 . (14) 

For the CKM matrix elements, we adopt the Wolfenstein parametrization and the updated pa¬ 
rameters A = 0.814, A = 0.22537, p = 0.117 ± 0.021, and p = 0.353 ± 0.013 [1]. 

Utilizing the above chosen distribution amplitudes and the relevant input parameters, we can 
get the numerical results in the pQCD approach for the form factor ^ at maximal recoil as 

follows, 

= 0) = . ds) 

where the errors arise from the shape parameter ojh in B meson distribution amplitude, the scalar 
decay constant fa^, and the Gegenbauer moments B^°{i = 1, 3) in the light oq distribution ampli¬ 
tude, respectively. This value agrees well with O.dO^Q os as given in Ref. [27]. The tiny deviation 
is just from the zero vector decay constant fao assumed in this work. 

B. CP-averaged branching ratios and CP-violating asymmetries 

In this subsection, we will analyze the CP-averaged B a^ao branching ratios and the CP- 
violating asymmetries in the pQCD approach at LO level. For B —)■ aoOo decays, the decay rate 
can be written as 

r = ^||^(l-2r2 ,)|d(B^o„o„)p, (16) 

where the decay amplitudes A can be referred correspondingly in Eqs. (6-8). Using the decay 
amplitudes obtained in last section, it is straightforward to numerically evaluate the CP-averaged 
branching ratios with errors as collected in Eqs. (17)-(19), 

Bt{B^ ^ ao+ao") = 1.5;0;^(a;,);°J(/,J10;e^(i?“o)+o.i(CKM) x , (17) 

Br{B+ ^ ao+ao°) = 6.li;?(a;fe);i;^(AJii(i?“°);n(CKM) x lO"® , (18) 

Bt{B^ ^ ao°ao°) = 2.7;j;J(a;,);°;®(Aj;j:3(i?“o);oi(CKM) x ; (19) 

^ The form factor can be extracted directly from Eq. (29) in [9] with the state S being ao. Of course, the 

readers can also refer to Ref. [27] for more details. 
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The dominant errors are indueed by the uncertainties of the shape parameter Ub = 0.40 ± 0.04 GeV 
for B meson, the scalar decay constant and the Gegenbauer moments B1°{i = 1,3) for the 
scalar ao(see Eq. (12) for detail), respectively. It is worth stressing that the effective constraints on 
the above mentioned non-perturbative parameters might be helpful to explore the QCD dynamics 
involved in these decays and to reveal the inner structure of the light scalar oq state. 

From Eqs. (17)-(19), one can obviously observe that the large B —)■ aoao decay rates are in 
the order of 10“® ~ 10“^ calculated in the pQCD approach at EO level, which could be easily 
detected through the dominant oq to r]7r(or tttt) final state [28] at the running EHC and the forth¬ 
coming Super-B/Belle-II experiments. As mentioned in the Introduction, some decays involving 
scalar mesons were suggested as the ideal channels to test the validation of the factorization as¬ 
sumption [14]. It is therefore worth stressing that the B'^ —)■ Og Og mode would be the best choice, 
because it only contains a significantly suppressed factorizable emission contribution and a neg¬ 
ligible non-factorizable emission contribution as proposed in naive factorization, but has a large 
branching ratio that could be easily tested in the near future experiments. Therefore, the obser¬ 
vation of this large B^ —)■ decay rate, on one hand, could offer an effective test to the 

breaking effects of the factorization assumption; on the other hand, might verify the qq compo¬ 
nents of the light scalar og evidently. Furthermore, it is surprising to find that the conventionally 
so-called ’’color-suppressed” B^ —OgOg mode has the largest branching ratio as 2.7 x 10“^, 
which is highly different from the known color-suppressed B PP modes, such as the famous 
B^ —)■ 7r°7r° channel with very small branching ratio around (9(10“^), although they embrace 
the same components at the quark level. Consequently, the hierarchy of the branching ratios ex¬ 
hibits theoretically as Br{B^ Oq^o) ~ Br{B^ —)■ Og Og ) > Br{B^ —)■ Og Og) in the pQCD 
approach, which is also dramatically different from that in the B ^ ttit decays as Br{B^ — 
^ Br{B~^ —)■ 7r+7r°) >> Br{B^ 7r°7r°) within theoretical errors [6, 12, 24, 25] and 
Br{B^ —)■ 7r+7r°) > Br(B^ —tt+tt”) > Br{B^ —tt^tt^) within experimental uncertain¬ 
ties [1, 2], respectively. In terms of the central values of the B —aoao decay rates, the following 
relation can be easily found, 

Br{B^ a°a°) > Br{B^ a+Ug ) > Br{B+ a+a°) , (20) 

which can be traced back to the factorization formulas as given in Eqs. (6)-(8). Specifically, the 
tree dominant contributions of these three decays are C 2 {Mnfa — Mnfs), Ci Mnjs + C 2 Mnja, and 
(Cl + C 2 )Mnfs, respectively, in which C 2 is much larger than Ci in magnitude with C 2 ~ 1.12 
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and Cl ~ —0.27 at the scale, and Mnfs{Mnfa) stands for the amplitude of the non-factorizable 
emission (annihilation) diagrams induced by the tree operators Oi 2 - The underlying reason is 
that, as presented in Eq. (11), the asymmetric leading twist distribution amplitude turns 

the originally destructive interferences induced by the symmetric one between the two 

non-factorizable emission diagrams, namely. Fig. 1(c) and 1(d), in the B —)■ PP decays into the 
presently constructive ones in the B ^ aoao modes. Meanwhile, the analogous phenomenon also 
occurs in the annihilation topologies. Note that the values of Mnja are usually a bit smaller than 
those of Mnfs in modulus, because the former is always power l/ms suppressed with ms being 
the B meson mass. It is interesting to note that the QCD behavior in light scalar oq is greatly 
different from that in the pseudoscalar pion, which can be seen apparently that the leading twist 
ao(pion) distribution amplitude is governed by the odd(even) Gegenbauer polynomials [10,29, 30]. 
Therefore, large non-factorizable contributions are observed in the 5 —)■ aoOo decays. 

In view of the surprisingly large Br(B^ —)■ a^al) and the amazingly small Br{B^ —)■ 7r°7r°) 
in the pQCD approach at LO level, respectively, we here present the numerical decay ampli- 
tudes^(See Tables I and II for detail) arising from every topology to clarify the aforementioned 
predictions explicitly. It can be clearly seen that the decay amplitudes in the B —)■ aoao decays 
exhibit very different pattern from those in the 5 —)■ tttt ones, although they embrace the same 
diagrams at the quark level: the former modes determined by the non-factorizable contributions 
with a larger scalar decay constant fa^ ~ 0.365 GeV, while the latter ones dominated by the factor- 

TABLE I. The factorization decay amplitudes(in unit of 10“^ GeV^) of the charmless hadronic B —)• 
aoao decays in the pQCD approach at leading order level, where only the central values are quoted for 
clarification. 


Decay modes 

Afs 

'^nfs 

■^nfa 

•^fa 

a([ao 

0.950 - f0.390 

1.619 - i2.982 

-1.056 - il.876 

-0.044-h fl.212 

B+ aj-af] 

-0.018 -h iO.007 

-1.268-hf2.926 

0.0 

0.0 

^ agag 

0.691 - to.284 

2.458 - i5.100 

-0.799 - il.363 

-0.035 -7 fO.853 


^ The topological amplitudes A/s, Anfs,Anfa, and Afa shown in the Tables I and 11 stand for the decay amplitudes 
of factorizable emission, non-factorizable emission, non-factorizable annihilation, and factorizable annihilation 


diagrams, respectively. 
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TABLE II. Same as Table I but for the charmless hadronic B ^ mr decays. 


Channels 

Afs 

'^nfs 

'^nfa 

Afa 

1 

+ 

t 

O 

cq 

-1.845 - i2.957 

0.095 -h i0.075 

-0.047-h iO. 159 

0.038-hiO.196 

B~^ —>■ 7r'*'7r^ 

-0.844- il.988 

-0.086 - f0.082 

0.0 

0.0 

fiO ^ 7r°7r° 

-0.461 - i0.104 

0.153-fiO.135 

-0.033-hiO.113 

0.029-hiO.139 


izable emission contributions with a smaller ~ 0.130 GeV, apart from the special —)■ 

channel. As mentioned above, the underlying reason is that these considered modes include dra¬ 
matically different QCD dynamics. Notice that, for the B —)■ aoOo decays, because of the vanished 
vector decay constant fa^ ~ 0, Ajs come only from the penguin contributions induced by the 
{S + P){S — P) operators, which are from the {V + A){y — A) ones by Fierz transformation. 
However, the phenomenologies shown in i? ^ aoao decays indicate that the famous B ^ tttt 
puzzle could be resolved if a new QCD mechanism is resorted to enhance the non-factorizable 
contributions. Of course, it is nontrivial to resolve the i? —)■ tttt puzzle just by including the large 
non-factorizable contributions. This point has been clarified in the literatures, for example, see 
Refs. [12,31]. 

Because of the large errors induced by the much less constrained hadronic parameters such as 
the scalar decay constant fa ^, the Gegenbauer moments Bi and B^ in the Oq distribution ampli¬ 
tudes, we derive the ratios of the branching ratios, in which the parameter uncertainties may be 
greatly canceled and be more helpful for measurements in the relevant experiments. 


Ro+ 

Roo 

R+0 


Br{B^ — )■ Oq Oq ) 
fir(5 + —)■ Oq Oq) 
Br{B^ —)■ Oq Oq ) 
fir(fiO ^ ajjag) 
Br{B~^ —)■ Oq Oq) 

Br{B0 a[]a[]) 


^ 0.56inj(u;,)in°(/,Jinj(fir)in°(CKM) , 

^ 0.23lHo(^.)^Ho(/ao)^Ho(i?r)^Ho(CKM) ; 


( 21 ) 

( 22 ) 

(23) 


It is well known that the fi —)■ tttt modes can provide important information to constrain the 
CKM unitary angle a. As they contain the same quark diagrams as the fi —)■ tttt decays, it is 
generally believed that the fi aoao processes can also provide complementary constraints on 
the angle a. Here, we show the a dependent branching ratios of the fi —?■ aoao decays in the 
pQCD approach at the LO level. Based on Eqs. (6)-(8), the decay amplitudes of fi —)■ aoao decays 
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can be rewritten as follows, 

^ = V:,VudT - V*VtdP = , 


(24) 


where the weak phase a = arg 


VtlVtd 

V*)^Vud 


, the ratio z = \V^lVtd/V*f,Vud\ ■ \P/T\, and 5 is the 


relative strong phase between tree(T) and penguin(P) amplitudes. Correspondingly, the decay 
amplitudes of the B aoao decays can be read as. 


A = KiViT - V,tV;,P = KiV;;T(1 + ze-'-^+‘>) , 

Therefore, the CP-averaged branching ratio of the B —)■ aoao decays shall be the following. 


(25) 


Br(B ->■ aoao) = (|^|^ + \ A\^)/2 = |lC),14dT|^(l + 2z cos a cos 5 + z^) . 


(26) 


It is thus easy to see that the CP-averaged branching ratio is a function of cos a for the given ratio 
z and the strong phase 5, which can be perturbatively calculated in the pQCD approach. This gives 
a potential method to determine the CKM angle a by measuring the CP-averaged branching ratios 
with precision. The dependence on the CKM weak phase a of the CP-averaged branching ratios 
for 5° —)■ OoUq (Solid line), —)■ Og ao(Dashed line), and B^ —)■ aoao(Dash-dotted line) decays, 

respectively, are presented in Fig. 2, where the central values of the predictions in the pQCD 
approach are simply quoted for clarification. Then we can directly observe that the central decay 
rates for the B —)■ agOo decays in the pQCD approach at LO level correspond to the value around 
90° of the CKM angle a, which agrees well with the constraints from various experiments [1]. 



FIG. 2. (Color online) Dependence on the CKM angle a of the i? —o^Og (Solid line), a^ag(Dashed line), 
and agOglDash-dotted line) decay rates at leading order in the pQCD approach, respectively. 

Now we turn to the evaluations of the CP-violating asymmetries of i? ^ aoao decays in the 
pQCD approach. For —)■ Og Og decay, the direct CP-violating asymmetry Aqp can be defined 
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as: 




( 27 ) 


Using Eq. (27), it is easy to calculate the direct CP-violating asymmetry for the considered B'^ —)■ 
Oq Oq mode as listed in Eq. (28), 

^ oJfflS) = -0.6iS:5(w,)lS:5(B“")lS:;(CKM)%, (28) 

This tiny direct CP-violating asymmetry would be hard to be measured because of the extremely 
small penguin contributions in magnitude, although the large strong phase can be obtained due 
to the constructive interferences between the two non-factorizable emission diagrams with the 
asymmetric oq leading twist distribution amplitude, which is very different from that in the —)■ 
7r+7r° mode with the small non-factorizable emission contributions, relative to the purely real 
amplitudes from the factorizable emission diagrams in the pQCD approach at EO level. 

As to the CP-violating asymmetries for the neutral decays B^ aoo-o, the effects of ^ 
mixing should be considered. The CP-violating asymmetries of B^(B^) —)■ ao’''ao“ and ao°ao° 
decays are time dependent and can be defined as 

r (b (At) /cp) — r {B^{At) —)■ /cp) 


Acp = 


r [B^At) ^ /cpj + r {B%At) ^ /cp) 

= cos(AmAf) -f A™p''sin(AmAf), 


(29) 


where Am is the mass difference between the two B^ mass eigenstates. At = tcp — Uag is 
the time difference between the tagged B^ (B^) and the accompanying B^ (5°) with opposite b 
flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct- and mixing-induced 
CP-violating asymmetries A^p and A™p^ can be written as 


4dir _ I-^CpI^ - 1 

i + iAcpr 


^mix _ 

/icp — 


2Im(Acp) 

1 + |AcpP’ 


with the CP-violating parameter A 


CP 


X VtlVu {fcp\Heff\B) 

^CP = Vf 


(30) 


(31) 


where ? 7 / is the CP-eigenvalue of the final states. Then the direct- and mixing-induced CP-violating 
asymmetries for the 5° —)■ and aJ^aQ decays in the pQCD approach at LO level can be 

calculated as, 

^ a+a„-) = , (32) 

^ a+o„-) = 0.9«:l(u,,)iJ‘(Br)ls:6(CKM)% , (33) 
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^ «o«o) = 16.2;J:I(a;,);i:®(i?r)lSi(CKM)% , (34) 

^ ay,) = 4.6l^:^(a;,)ll|(i?“o)+|9^cKM)% , (35) 


where we have negleeted the vanishing theoretieal errors for the CP-violations in i? —)■ aoOo de¬ 
cays arising from the scalar decay constant fa^ of oq meson. It is interesting to see that these 
two channels, namely, —)■ and B^ —)■ a^a,, have large branching ratios and large direct 

CP asymmetries simultaneously, which could be easier to be measured at the running LHC ex¬ 
periments and the forthcoming Super-B/Belle-II factory, and have the potential to reveal the QCD 
dynamics and the inner structure involved in the light scalar oq meson. 

Similarly, based on Eqs. (24), (25), and (28), the direct CP-violating asymmetry can also be 
expressed as the function of the CKM angle a. 


2z sin a sin S 
1 + 2z cos a cos 5 + z"^ 


(36) 


Then the precise measurements on these large direct CP violations can also provide the constraints 
on the CKM angle a potentially. The variation of the direct CP-violating asymmetries with the 
CKM angle a for the 5° —)■ ag'ao (Solid line) and aoao(Dashed line) decays is shown in Fig. 3. 
Again, the central value about 90° of the CKM angle a can be utilized to produce the above 
mentioned large direct CP violations. 



FIG. 3. (Color online) Dependence on the CKM angle a of the B ^ a, a, (Solid line) and agaglDashed 
line) direct CP violations at leading order in the pQCD approach, respectively. 


IV. SUMMARY 

In summary, we studied the two-body charmless hadronic B —)■ aoao decays, which have the 
same Feynman diagrams as the i? —)■ tttt modes at the quark level, by employing the pQCD factor- 
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ization approach based on the faetorization theorem. Based on the assumption of two-quark(gg) 
strueture of the light sealar oq state, we make theoretieal predietions on the CP-averaged braneh- 
ing ratios and the CP-violating asymmetries of the eonsidered B —)■ aoao ehannels in the SM. Due 
to the large non-faetorizable eontributions indueed by the asymmetrie leading twist distribution 
amplitude of oq meson, large branehing ratios in the order of 10“® ~ 10“^ have been predieted 
in the pQCD approaeh at LO level. At the same time, large direet CP violations around 15% and 
30% in the —)■ and deeays have also been observed. It is therefore expeeted that 

the large branehing ratios plus the large CP asymmetries would be easier to be measured at the 
running LHC experiments and the fortheoming Super-B/Belle-II faetory, if oq is indeed the qq 
bound state. Furthermore, the large non-faetorizable eontributions in the B —)■ aoao deeays ean 
hint some important information on resolving the famous i? —)■ tttt puzzle, although this is non¬ 
trivial work as elarified in the literatures [12, 31]. The deteetion of these eonsidered deeays might 
be helpful to investigate the QCD dynamies in the ehannels and to explore the inner strueture of 
the light sealar oq state. The investigation of the B —)■ aoao deeays eould also provide more eom- 
plementary eonstraints on the CKM weak phase a, sinee the same eomponents as the i? —)■ tttt 
modes exist in the eonsidered B aoao ones at the quark level. Frankly speaking, the predie¬ 
tions in the present work suffered from large uneertainties indueed by the mueh less eonstrained 
hadronie parameters sueh as the Gegenbauer moments B^° and B^°, whieh need further studies in 
the non-perturbative QCD(sueh as QCD sum rule and/or Lattiee QCD) ealeulations and the rele¬ 
vant experimental measurements(e.g., at BESIII, LHC, Super-B/Belle-II, ete.) on the produetions 
and/or deeays involving the oq state. 
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